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Abstract 

The dynamics of generic stochastic Lotka-Volterra (discrete logistic) systems 
of the form Q] Wi(t + 1) = X(t)wi(t) + aw{t) — bwi(t)w(t) is studied by com- 
puter simulations. The variables Wi, i = 1,...N, are the individual system 
components and w(t) = jfJ2i w i(t) 1S their average. The parameters a and b 
are constants, while \(t) is randomly chosen at each time step from a given 
distribution. Models of this type describe the temporal evolution of a large 
variety of systems such as stock markets and city populations. These systems 
are characterized by a large number of interacting objects and the dynamics is 
dominated by multiplicative processes. The instantaneous probability distri- 
bution P(w, t) of the system components Wi, turns out to fulfill a (truncated) 
Pareto power-law P(w, t) ~ w^ 1 ^'*. The time evolution of w(t) presents inter- 
mittent fluctuations parametrized by a truncated Levy distribution of index 
a, showing a connection between the distribution of the w^s at a given time 
and the temporal fluctuations of their average. 
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I. INTRODUCTION 



Power-law distributions have been observed in all domains of the natural sciences as well 
as in economics, linguistics and many other fields. Widely studied examples of power law 
distributions include the energy distribution between scales in turbulence [IJ , distribution of 
earthquake magnitudes 0, diameter distribution of craters and asteroids ||, the distribution 
of city populations 0,0, the distributions of income and of wealth |HTI|, the size-distribution 



of business firms [n,l2"| and the distribution of the frequency of appearance of words in texts 
H. A related phenomenon is the fact that in a variety of systems the temporal fluctuations 
exhibit a scale invariant behavior in the form of (truncated) Levy-stable distributions. Well 
known examples are the fluctuations in stock markets |7|,|T^]. 

Although systems which exhibit power-law distributions have been studied extensively in 
recent years there is no universally accepted framework which can explain the origin of the 
abundance and diversity of power-law distributions. One context in which the emergence of 
scaling laws and long range correlations in space and time is well understood is equilibrium 
statistical physics at the critical point |14| - |r^] . By contrast, scaling behavior, power law 
distributions as well as spatial and temporal power law correlations in generic natural 
systems is still the subject of intense study |]T8|-|31|1 . 

An approach that proved to be useful in the study of complex systems is to identify 
for each system the relevant elementary degrees of freedom and their interactions and to 
follow-up (by monitoring their computer simulation) the emergence in the system of the 
macroscopic collective phenomena [^] . This approach was applied to the study of multiscale 
dynamics in spin glasses |33] and stock market dynamics f34|. Using a generic class of 
models with a large number of interacting degrees of freedom, it was shown that macroscopic 
dynamics emerges under rather general conditions. This dynamics exhibits power law scaling 
as well as intermittency [j34H36|| . These models turn out to be particularly suitable to describe 
systems such as stock market dynamics with many individual investors [p?|- |39|j4l|j4^j42| 



where each system component describes a single investor (or stock pi). Such systems 
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involve complex temporal dynamics of many degrees of freedom but no spatial structure. The 
models introduced in ]36| can also describe systems such as population dynamics |P6|- ^9| , ^3| fl , 



spatial domains in magnetic or turbulence models [pT] , p2| or regions in generic phase 
spaces p3^p5[, which have spatial dependence. 



In this paper we present numerical studies of generic stochastic Lotka-Volterra systems. 
These systems |36| basically consist of coupled dynamic equations which describe the discrete 



time evolution of the basic system components Wi, i = 1, . . . , N. The structure of these 
equations resembles the logistic map and they are coupled through the average value w(t). 
The dynamics includes autocatalysis both at the individual level and at the community 
level as well as a saturation term. We find that under very general conditions, the system 
components spontaneously evolve into a power-law distribution P(w,t) ~ w~ l ~ a . The 
time evolution of w(t) presents intermittent fluctuations parametrized by a (truncated) 
Levy-stable distribution with the same index a, showing an intricate relation between the 
instantaneous distribution of the system components and the temporal fluctuations of their 
average. 

The paper is organized as follows. In Sec. II we present the generalized logistic model 
introduced in P5[ . Simulations and results are reported in Sec. III. Discussion of previous 



results as well as of our findings is given in Sec. IV, and a summary in Sec. V. 



II. THE MODEL 



The generalized logistic system |36[ describes the evolution in discrete time of iV dynamic 
variables w^, i = 1, . . . , N. At each time step t, an integer i is chosen randomly in the range 
1 < i < N, which is the index of the dynamic variable W{ to be updated at that time step. 
A random multiplicative factor A(t) is then drawn from a given distribution 11(A), which 
is independent of i and t. This can be, for example, a uniform distribution in the range 
A m in < A < X max , where \ m in and \ ma x are predefined limits. The system is then updated 
according to 



3 



Wi (t + 1) = X(t)wi(t) + aw(t) - bwi(t)w(t) 

Wj (t + l) = Wj (t), j = l,...,N; j^i. (1) 

This is an asynchronous update mechanism. The average value of the system components 
at time t, is given by 

N 

N 



= 4 £«>*(*)• ( 2 ) 

%=i 



In general, using instead of the average, a weighted average of the w^s would lead to similar 
results. The parameters a and b may, in general, be slowly varying functions of time, however 
we will now consider them as constants. The first term on the right hand side of Eq. (0), 
describes the effect of auto-catalysis at the individual level. For instance, in a stock-market 
system it represents the increase (or decrease) by a random factor X(t) of the capital of the 
investor % between time t and t + 1. The second term in Eq. (|l|), describes the effect of 
auto-catalysis at the community level. In an economic model, this term can be related to the 
social security policy or to general publicly funded services which every individual receives. 
In molecular or magnetic systems, this term may represent the mean-field approximation to 
the effect of diffusion or convection ||50|| . The third term in Eq. (|I|), describes saturation 
or the competition for limited resources. In an ecological model, this term implies that 
for large enough densities, the population starts to exhaust the available resources and 
each sub-population loses from the competition over resources a term proportional to the 
product between the average density population and its own size. We refer to Eq. ([!]) as the 
generalized discrete logistic (GL) system because when averaged over i, this system gives 



the well known discrete logistic (Lotka-Volterra) equation |56|]5? 



w (t + 1) = (X + a)w(t) -bw 2 {t). (3) 

In the general case, the parameters a, b and the distribution 11(A) may depend on time. 
Consequently, even the solution of the asymptotic stationarity condition w(t + 1) = w(t) 
may depend on time according to 
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w{t) = {\{t) + a-l)/b{t). (4) 

In fact, the typical dynamics of microscopic market models [^-^T| is generically not in a 
steady state. As will be shown below, systems which exhibit an effective GL dynamics [Eq. 



(1)] lead f36[, under very general conditions, to a power law distribution of the values wf 

P(w) ~ w- l ~ a . (5) 

Moreover, the time evolution of w(t) presents intermittent fluctuations following a (trun- 
cated) Levy-stable distribution with the same index a, 



III. SIMULATIONS AND RESULTS 

To examine the behavior of the GL model introduced in ||36j| , Eq. ([!]) we performed 
extensive computer simulations. Most simulations were done with N = 1000 system compo- 
nents, using various values of the parameters a and b and different distributions 11(A) of the 
multiplicative factor A. We focused on the power law distribution of the system components 
Wi as well as on the fluctuations of w. Fig. [I] shows the distribution of Wi, i = 1, . . . , N, 
obtained for iV = 1000, a = 0.00023, b = 0.01 and A uniformly distributed in the range 
1.0 < A < 1.1. A power law distribution is found within the range 

w < Wi < Nw (6) 

which is bounded from below by the average wealth and from above by the total wealth, 
and spans nearly three decades. 

The robust nature of the power-law distribution is demonstrated in Fig. ^ for 6 = 0. In 
this case w(t) does not reach a steady state and keeps increasing (or decreasing) indefinitely. 
However, the power-law behavior is maintained. Moreover we find that the the exponent 
a is insensitive to variations in b: even for values of b differing by an order of magnitude 
(corresponding to w varying by an order of magnitude), the power law exponent a is virtually 
unchanged. 
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At each time step the system component to be updated is chosen randomly. Since the 
system components Wi exhibit a power- law distribution, given by Eq. (|5|), the impact of the 
update move on w{t) exhibits a broad distribution. The dynamics involves, according Eq. 
(|I|), a generalized random walk with steps sizes distributed according to Eq. (|5]). Therefore, 
the stochastic fluctuations 

w{t + r) - w{t) 
{ ' w{t) V) 

of w(t) after r time steps, are governed by a truncated Levy stable distribution L a (r). 
This means that rather than shrinking like iV" 1 / 2 the fluctuations of w(t) have infinite 
variance in the thermodynamic limit (modulo the truncation). The truncation in the Levy- 
stable distribution corresponds to the cutoff's in the power-law distribution, given in Eq. 

Typically, the truncation in r is bounded by the relative width of A times the largest 
Wi/{Nw) value. 

Fig. shows the distribution P(r) of the stochastic fluctuations r(r), for r = 50, which is 
given by a Levy-stable distribution L a {r). We find indeed that all the values of r are smaller 
than the relative width of A (0.1) times the maximal value of Wi/(Nw) from Fig. [I]. The 
cut-off in the distribution of the temporal fluctuations originates therefore in the cut-off in 
the Pareto power law in Fig. |l|. In the absence of this cut-off the variance of the distribution 
of fluctuations would be infinite. The divergence of the variance modulo finite size effects is 
analogous to the divergence of the susceptibility in ordinary statistical mechanics systems 
at criticality. 

The peak of the (truncated) Levy-stable distribution scales with r according to 

L a (r = 0) - t- 1 ^ (8) 

where a is the index of the distribution. In Fig. || we show the height of the peak P(r = 0) 
of the distribution of fluctuations in w as a function of r for the parameters used in Fig. [I], 
which give rise to a power-law distribution of the w^s with a = 1.4. It is found that the slope 
of the fit in Fig. |] is —0.71 which is equal to following the scaling relation of Eq. (||). 
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This is a further indication that the fluctuations of w follow a Levy-stable distribution with 
the index a which equals the exponent of the Pareto power law in Fig. |1]. It is gratifying 
that an explanation of the 100 years old Pareto power law in these non-equilibrium systems 
which we have studied is provided by a straightforward explicitation of the almost as old 
Lotka- Volterra equation [0,0 . 

To provide more intuition about the dynamics leading to the power law distribution 
of Wi, we show in Fig. (|5|) the time evolution of a GL system starting from a uniform 
distribution of Wi, i — 1, . . . , 1000. We observe that the distribution gradually broadens. In 
the first stages it becomes of log-normal form and then it evolves into a power-law as the 
non- multiplicative effects in the vicinity of the lower bound become significant. 



IV. DISCUSSION 



A. Previous Results 



The numerical results of the previous Section show convincingly that generic (even non- 
stationary) systems with effective dynamics governed by the GL system of Eq.(|l|), lead to 
(truncated) Pareto distributions of the system components. They also lead to (truncated) 
Levy-stable laws of the fluctuations of the average. Let us now explain intuitively why this 
is the case. Consider first the Kesten system which is well known to present power laws 



w{t + 1) = \{t)w{t) + p{t) (9) 

where the random numbers A and p are extracted from two positive distributions independent 
of t. The Kesten system has a number of shortcomings which makes it unfit for most practical 
applications in natural systems: 



In Eq. (^|), there is only one variable, (no index i). It describes a non- interacting 
investor (animal, city) in a market (ecology, country) which induces effectively to him 
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the return (growth) \(t) — 1 after each trade (reproduction/replication/multiplication) 
period t. 

• In order for this system to exhibit a power law distribution, A has to be predominantly 
less than 1 such that it causes Xw to be on average smaller than w. Otherwise, the 
resulting w distribution is a log-normal with width expanding in time. This would 
correspond in the infinite time limit to a power law of the form P(w) ~ w~ l . The 
dependence of the Kesten model on a shrinking dynamics is incompatible with most 
of the natural systems in which the growth is positive. For instance, the shrinking 
multiplicative dynamics is certainly not a good model for a stock market where the 
investors i expect their wealth to increase on average (otherwise they just stay out of 
the market). 

• In realistic markets (ecologies, societies), the average wealth (population) w varies sig- 
nificantly in time. In the Kesten model this can be realized only by varying the 
distributions of A and p which in turn would significantly affect the exponent a of the 
power law [Eq. (||)]. In the GL system, on the other hand, changes in the environment 
are represented by changes in the coefficient b of the resources limitation/competition 
term. This can lead to changes by orders of magnitude in the total wealth/population 
Nw without affecting the exponent a. Interestingly, it turns out that the exponent a, 
in the distribution of wealth, has been stable for the last one hundred years and across 
most western (capitalist) countries. 

We will see later how the GL system solves the shortcomings of the Kesten model. 
Meantime let us give an intuitive explanation of why the Kesten system leads to a power 
law given by Eq. (|5|). First one should realize that due to the p(t) term in Eq. (|9|), the 
values of w(t) are typically kept above a certain minimal value of order p. Let us therefore 
effectively substitute the p term in the Kesten equation [Eq.(|9])] with the condition that 
w(t) > p. More precisely, each time w(t) becomes smaller than p, it is reposed "by hand" 
to the value p. 
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In the resulting system: w(t + 1) = X(t)w(t), with w(t) > p one can take the logarithm: 
lnw(t + 1) = lnw(t) + ln\(t). The lower bound condition becomes then lnw(t) > Inp. This 
represents a system in which Inw undergoes a random walk with a drift towards smaller 
values and with a reflecting barrier at Inp. One can compare this with a molecule in grav- 
itational field submitted to the collisions with the rest of the gas (resulting in friction and 
Brownian motion) and bounded from below by the ground level. It is not surprising there- 
fore that (by analogy to the barometric equation) the resulting probability distribution for 
Inw is an exponential: 

p(lnw) ~ e~ pinw (10) 

which written in terms of w itself gives: 

P(w) ~ w' 1 ' (11) 

The particular value of (3 depends in the Kesten system on the details of the distribution 
11(A) and is such that the drift towards lower values induced by A is balanced by the drift to 
larger values induced by p. As a consequence, this model, if (mistakenly) applied to the stock 
market (ecology, society etc.), would predict not only negative average returns (growth) but 
also an exponent a in the power law that is highly sensitive to the parameters [p2 H64lj6q,|67| . 
On top of all these shortcomings, the Kesten system does not predict a (truncated) Levy- 
stable distribution of the w fluctuations (as repeatedly measured in nature |19|-p^]). To get 
the (truncated) Levy-stable distribution the following conditions should be satisfied: the 
index i, of the component W{ to be updated at time t, is chosen randomly, the w^s satisfy 
a power-law distribution and the update step is multiplicative, namely, the change in Wi is 
proportional to its current value. For example, even if the dynamics leads to a power-law 
distribution of Wi, the fluctuation may not be described by Levy-stable distribution if the 
magnitude of the update of Wi is not proportional to u>, itself. 
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B. How Does Our Model Work 



Let us now see how the GL model [36| solves the problems with the Kesten system. The 
main new ingredients in the GL model are the appearance of w and the appearance of an in- 
dex i in u>j. These two objects allow the introduction of a crucial ingredient which was absent 
in the Kesten system: the interaction between the investors (sub-ecologies, sub-systems). 
While the interaction in the stock market (ecology, society) is represented in the Kesten 
system only implicitly by the stock returns (growth) X(t) — 1 we introduce now additional 
interactions between the investors (individuals, families) i which are mediated by the aver- 
age w(t) and are crucial for the dynamics of the system. Obviously, such terms, containing 
w(t) could not appear in an equation like Kesten which considers only the dynamics of one 
variable at a time. In order to introduce the crucial terms including w one has to give up 
the picture of a single random investor and to embrace the picture of a macroscopic set of 
microscopic investors, interacting among themselves through the market mechanisms. The 
result is the system of nonlinear GL equations which are coupled through w(t) [Eq. (p])]. In 
order to gain insight into the emergence of the power law and Levy-stable intermittency in 
the GL system, one can express it formally as: 

Wi (t + 1) = [A(t) - b{t)w{t)} Wiit) + a{t)w{t) (12) 

If one ignores for the moment the effect of the changes of w^s on the value of w, the system 
[Eq. ([T2|) 1 is of the Kesten type [Eq. (H)] and we expect therefore the emergence of a scaling 
law, given by Eq. @. If the effect of the changes in Wi on w is considered, then one sees 
that (for non- vanishing b) the system is self-tuning towards the value of w given by Eq. ([|). 
This self-tuning is realized by the dynamics of the average in Eq. (0). If w(t) is small, then 
according the first term in Eq. fll2|), Wi will typically increase and will make w(t) increase 
too. If w(t) is large, then according the first term in Eq. flT2|), Wi will typically decrease and 
will make w{t) decrease. While in the synchronous Lotka-Volterra [with a global time step 
updating based on Eq. (3)] the system may have large steps and get into behavior alternating 
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chaotically between large and small w(t) values, in the case of the sequential updating [Eq. 
(1)] of the WiS, the average will eventually self tune to a value of w(t) given by Eq. (|]). 
The fluctuations around this value will be dominated by the first term in Eq. (|T2|) and will 
consist of a random walk with steps proportional to u>j. Since W{ are distributed by a power 
law, the fluctuations will be distributed by a Levy-stable distribution of corresponding index 
|69| , [70|| . In order to understand why the W{ distribution is only weakly dependent on b(t), 



one can substitute Eq. (§]) in to Eq. fll^) and use the normalized variables V{ = Wi/w. One 
then obtains an equation of the Kesten form: 

Vi(t + 1) - Vi(t) = \X(t) - X(t) - a(t)l Viit) + a(t) (13) 

Note that we used here the approximation that the dynamics of w is much slower than the 
dynamics of W{. One observes that both b(t) and w(t) are absent from Eq. (pT3[): their 
respective effects cancel. In fact, one finds in simulations (Fig. 0) that the distribution of 
Vi(t) (and therefore of Wi(t)) fulfills a power law of Eq. (|5|) with exponent a independent 
on the variations of w(t). One also sees from Eq. ( |T3"D that the dynamics is invariant to an 
overall shift in the distribution 11(A). This means that in particular the GL multiplicative 
factors A can be significantly (and generically) larger than unity allowing (in contrast to the 
Kesten system) for expanding (growing) dynamics. Eq. (|l^) implies time correlations in the 
amplitude of the fluctuations of w. It was brought to our attention by D. Sornette that our 
data seem consistent with the log-periodic corrections due to complex exponents discussed 
in |)5| as well as other the other data collected in the stock market |T9|- p8| , [7T| . [7^ . [75| , |72| , [76 . 



Our mechanism relates the emergence of power-laws and macroscopic fluctuations to the 
existence of auto-catalyzing sub-sets in systems composed of many microscopic entities. In 
particular, the use of the w is not mandatory: generic systems of the type Wi(t + 1) = 
J2j KjWjit) — J2j,kbijkWj(t)wk(t) may also present similar properties. At a more conceptual 
level, the challenge is to identify, in an as wide as possible range of natural systems, the 
elementary objects i, the degrees of freedom Wi associated with them and the GL interactions 
explaining in each case the emergence of scaling and intermittency. 
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V. SUMMARY 



In summary, we have studied the dynamics of a generic class of stochastic Lotka-Volterra 
(discrete logistic) systems introduced in |36|] using computer simulations. These systems con- 
sist of a large number of interacting degrees of freedom Wi(t), i = 1, ...N, which are updated 
asynchronously. The time evolution of each system component is dominated by a stochastic 
individual autocatalytic dynamics, in addition to a global autocatalytic interaction medi- 
ated by the average w(t), and a saturation term. These models describe a large variety of 
systems such as stock markets and city populations. We find that the distribution P(w,t) 
of the system components Wi, fulfills a Pareto power-law P(w,t) ~ The average 

w(t) exhibits intermittent fluctuations following a (truncated) Levy-stable distribution with 
the same index a. This intricate relation between the distribution of system components 
and the temporal fluctuations resembles the behavior of a variety of empirical systems. For 
example, it provides a connection between the power-law distribution of wealth in society 
and the fluctuations in the stock market which follow a (truncated) Levy-stable distribution. 
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FIGURES 

FIG. 1. The distribution of wealth Wi, i = 1,...,N (the number of investors P(w) possessing 
wealth w) for N = 1000 investors obtained from a numerical integration of Eq.([l|) with parameters 
a = 0.00023, b = 0.01 and 11(A) uniformly distributed in the range 1.0 < A < 1.1. The distribution 
(presented here on a log — log scale) exhibits a knee on the left hand side and a broad tail of power 
law distribution on the right hand side. This power law behavior is described by P(w) ~ w^ 1 ^ ', 
where the exponent a = 1.4. The distribution is bounded by an upper cutoff around w max = Nw. 

FIG. 2. The distribution of the values of w h i = 1, . . . , N for N = 1000, a = 0.0001, b = 0.0 
and 11(A) uniformly distributed in the range 1.0 < A < 1.1. Because of the absence of the saturation 
term (6 = 0), the system is not stationary and w varies in time by orders of magnitude. In spite 
of this, the instantaneous normalized w distribution at each instant remains always a power law of 
constant exponent a. 

FIG. 3. The distribution of the variations of w after t steps r(r) = [w(t + t) — w(t)]/w(t), 
where r = 50, for the same parameters as in Figjl[ This distribution has a Levy-stable shape with 
a = 1.4. One can see that the shape on a semidogarithmic scale differs from a parabola (Gaussian 
distribution) in that it has significantly larger probabilities for large Wi values. 

FIG. 4. The scaling with r of the probability that r(r) = [w(t + t) — w(t)]/w(t) is 0. The 
parameters of the process are like in Fig. |l| and Fig. ||. The slope of the straight line on the 
logarithmic scale is 0.71 which corresponds to a Levy-stable process with a = 1/0.71 = 1.4. 

FIG. 5. The time evolution of P(w) for a system starting from an uniform distribution 
of Wi. In the first stages the distribution is log-normal and it then becomes power-law as the 
non-multiplicative effects at the lower bound start being effective. The process of convergence to 
the power law is much shorter than the actual equilibration of the w value. 
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